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I. INTRODUCTION 



We study the effective conductivity tensor of a two-dimensional composite consisting of 
a periodic array of circular cylinders embedded in a host matrix. The problem has been 
studied before by Rayleigh^ for the case of a square array of cylinders. His method was 
extended^ for regular arrays of circular cylinders. A method of functional equations^ 
employing analytic functions was used to find an expression of the conductivity tensor for 
small volume fraction of inclusions. For rectangular lattice of inclusions an efficient method 
based on the use of elliptic functions was suggested^, which is developed further in the 
present paper. Effective conductivity can be also evaluated numerically 7 . 

The goal of this work is to find an analytic expression for the effective conductivity tensor 
in the case of an arbitrary doubly periodic array of circular cylinders when the effective tensor 
can be anisotropic. The solution of the problem consists of two steps. First, we construct 
a quasiperiodic potential using a combination of Weierstrass ^-function and its derivatives. 
That ensures periodicity of the electric field in the whole plane and, as a result, avoids the 
problem of summation of a conditionally convergent series. This approach is similar to that 
in Ref. applied to biharmonic problems of the theory of elasticity. We reduce the problem 
to an infinite system of linear equations and find its solution in the form of a convergent 
power series (in terms of a parameter proportional to the volume fraction of the cylinders) 
whose coefficients are determined explicitly. Second, we determine the average electric field 
and the current density within one parallelogram of periods and find an exact expression of 
the effective conductivity tensor that relates the two quantities. 

II. DERIVATION OF PERIODIC POTENTIAL 

Suppose that a periodic lattice of identical circular inclusions of radius a is introduced 
into a uniform electric field E = [E x , E y ] T applied in the plane perpendicular to the cylinder 
axes. The nodes of the lattice in the complex plane are generated by a pair of vectors 
2(jj\ and 2u)2, Im ^ > (see Figured]). In polar coordinates the potential u(r, has the 
following properties: 



Au = 0, 




Uj n in the inclusion 



u ex in the medium 



(1) 
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FIG. 1. Circular inclusions of radius a arranged in a periodic lattice with periods 2co± and 2oj2- 



and on the boundary r = a 



(T; 



du in du ex 



(2) 
(3) 



'in r\ ^ ex n 1 

or or 

where <7j n and o ex are the electric conductivity of the inclusions and the medium, respectively. 
It is convenient to represent potential u in the complex form u = Re <&(z) with 



oo 

Ea^2(A n + iB n ) (- 



« 2n+l 




J 


\z\ 



z\ < a, 



n=0 



Vl a 2n+1 



$ ex (z) = (-E x + iE y ) z + EaJ2 ( C n + %D n ) C 



(2n), 



n=0 



z), \z\ > a, 



(4) 
(5) 



where A n , B n , C n , and D n are unknown real dimensionless coefficients, E = \E\, and ((z) 
is Weierstrass' ^-function 



1 1 : 

+ t; — + 



z — P ' P ' P 2 

<> 1 m,n 1 m,n 1 m,n 



(6) 



Here ^ 2n \z) denotes derivative of order 2n, and P m ,n — 2mo>i + 2na> 2 . Prime in the sum 
means that summation is extended over all pairs m, n except m = n = 0. 
Below we will use some properties of (^-function 9,10 , (^-function is an odd meromorphic 
function with simple poles at P mn . It has the quasiperiodicity property 



C(* + 2wi) = C(*)+27/i, 77i = C(wi), 
C(z + 2u 2 ) = C(z) + 2r] 2 , ri2 = C(u 2 ), 



(7) 
(8) 



where constants rji and 7] 2 are related by the Legendre identity 

Vi ^2 - V2 u\ = —. (9) 

Derivative of C,{z) is a periodic function and is expressed through Weierstrass elliptic function 

p(z) by 

C'(z) = -p(z). (10) 

This property ensures the electric field to be periodic in the medium, while (JT])-© guarantee 
that the potential changes by a constant value in the direction of either 2oj\ or 2uj 2 . 
Function p(z) satisfies the differential equation 

2 



dp(z) 



dz 

where g 2 and g 3 are two invariants defined by 



Ap 3 (z)-g 2 p(z)-g 3 , (11) 



^ = 60 E'75T-' 43 = 140^'-^, (12) 

n,m m ' n n,m m ' n 

which are used for numerical evaluation of p(z) and C( z )- 111 particular, we will use the 
following homogeneity property 

((z;g 2 ,g 3 ) = - e ({zr 1 ;g 2 £\g 3 f). (13) 

To satisfy conditions (E])- (JHJ) on the inclusion surface we expand ((z) and its even deriva- 
tives in a Laurent series 

C^W = ^-|^^^«-'. »>0, . = 0, (14) 

where 

a * = E / ^T' k = 2,3,.... (15) 
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Sums ( 1T51) contain only even powers of P m ,n since for every point P m ,n = 2mu\ + 2nu 2 on the 
lattice there exists symmetric point — P m .n and the sums with odd powers vanish. Also, if 
the periods 2ui\ and 2oj 2 of the lattice are fixed, sums Sk remain bounded as k — > oo. Thus, 
potential u(r, $) near the inclusion surface has the form 



Ea22(A n cos[{2n+ 1)0] - B n sin[(2n + 1)0]) (-) , (16) 

n=0 

/a \ 2n+l r 

-(^ cos t? + E y sin tf)r + £a ^- J [C„ cos(2n + + L> n sin(2n + 1)# 



n=0 



^ E io^f till * 2k+1 r 2n+1 \{C k cos(2n + l)d - D k sin(2n + 1)0) 



n+fc+l 



- (Cfc sin(2n + + D fc cos(2n + 1)0) si +k+1 , (17) 

where and s k denote the real and imaginary parts of the sum s k , respectively. 
From (J2])-© we obtain relations between the coefficients 

A n = ^— C n , (18) 

B n = D n , (19) 

and a system for determining C n and D n 



r — rv + 2m + 1)! , R I . 2n+2m+2 _ <■ (r)r >\ 

W a / J ^2m)\(2n + 1)! ' m n+m+1 Wm Vm+iJ " _ i? ,0 ' ^ ' 



m=0 



n ™ (2n + 2m + 1)! / r j n r \ 2n+2m+2 _ s / 91 \ 

D n -a (2m)\(2n+ 1)\ \ UmS ™+™+± ~ u ™ s n+m+i) a ~ ~E ,0 ' ^ ' 



() (2m)!(2n+l) 

with a = — — and 5 n fl being the Kronecker delta. 
Let us introduce the following notation: 



^ = min{2|wi|,2|w2|}, (22) 
h = j, (23) 



S * = Y!(jr-) = S k+ iS l, k = 2,3,..., Si = 0. 

n rn v ' ' 



(24) 



Here S^ and S k denote the real and imaginary part of S k , respectively. Note that once S2 
and S3 are computed, a recurrence relation*^ allows to evaluate higher order terms. The 
system then can be written as 
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(2n + 2m + l)\ , R _ 2 , 2n+2m+2 

"= / ^ (2m)!(2n+l)! ^ m "+ m+1 ^m&n+m+l) n 

y^ (2n + 2m + l)! f , _ fi U2n+ 

-t^n " / ^ (2m)\(2n + 1)\ °n+m+l ^m&n+m+l) " 



a — 5 nj0 , 



m=0 

or in vector form 



— a — o nfi , 



(25) 
(26) 



x n G n ^ m x m h 



2n+2m+2 



m=0 



where 



D r 



y 



a 
E 



En- 



E v 



G r 



(2n + 2m + 1)! 
(2n + l)!(2m)! 



"n+m+1 
' J n+m+l °n+m+l 



j ^0,0 — 0. 



(27) 



(28) 



(29) 



In Appendix we formulate conditions providing existence and uniqueness of the solution of 
the system ( 1271) in the space of bounded sequences, as well as the possibility to obtain its 
solution by truncation or in the form of a convergent power series in h. The latter method 
is used below. 

Let us look for solution of ( 1271) in the form of a power series in h 



2n+2m+2 



(30) 



m=0 



Substituting (|30|) into ( l2Tj) and equating the coefficients of like powers of h we obtain a 
recurrence relation for p n , m 



Pn,o = G nfi y, 

Pn,k ^ ^ G n ^ m Pm, 



k—2m—li 



fc = l,2,..., 



(31) 
(32) 



m=0 



where [s] denotes the integer part of s. Below we give several first terms of series expansion 
of x which are needed for calculation of the effective properties 



x = y + G ,xG 1>0 yh s + G , 2 G 2fi y h 12 + G Qjl G 1}1 G lfi y h 14 + ((G ,iGi, ) 2 + G , 3 G 3 ,o) y h 16 
+ (G ! o,iG ! i ) 2G ! 2 1 o + ^0,2^2,1(^1,0) y h 18 

+ (^0,1^1,0^0,2^2,0 + Go t iGi t iGx t Q + ^0,2^2,0^0,1^1,0 + (^0,4(^4,0) 2/ ^ 20 + 0{h 22 ) 

(33) 



or 



where 



a; = M ft y, 



(34) 



+ 30a J 



,11 



M & = (1 + 3a 2 |,S2| 2 /i 8 + 5a 2 \S 3 \ 2 h 12 ) I 

Si {[S 2 R ] 2 - [Si] 2 ) + 2S?SISZ Si ([Si] 2 - [Si] 2 ) - 2S*SIS£ 

Si ([Si] 2 - t^] 2 ) - 2^|^ -Si ([^] 2 - t^] 2 ) - 25?^ J 
+ (9a 4 |S 2 | 4 + 7a 2 |S 4 | 2 ) h 16 I 

qR I qR qR 1 ql ql\ 1 ql ( qR ql ql qR\ ql I qR qR 1 ql ql\ 1 C-R / ci?, 0/ 0/ cRN 
J 2 V 3 4 °3°4 / ' "r J 2 \°3 °4 °3°4 J J 2 \ J 3 °4 "r °3°4,/ ">" °2 \ J 3 °4 °3°4 ) 

ql I qR qR 1 q7 j_ ci? I qR ql ql qR\ ql I ql qR qR ql\ _ qR I qR qR 1 ql ql\ 
°2 \°3 °4 1* J 3 D 4) "1 D 2 \°3 °4 °3°4 ) °2 \ D 3 J 4: °3 °4) °2 \°3 °4 ° 3 4j 

S^ ([Si] 2 - [Si] 2 ) + 2S*SiSi Si ([S 2 R ] 2 - [Si?) - 2SiSiSi 
[ Si ([Si] 2 - [Si] 2 ) - 2SiSiSi -Si ([Si] 2 - [Si] 2 ) - 2SiSiSi 



+210a 



h 



IN 



+ 30a 6 



+ (30a 4 \S 2 \ 2 \S 3 \ 2 + 9a 2 \S 5 \ 2 ) l)h 20 + 0(h 22 ), 



(35) 



and / is the identity matrix. 

Matrix Mh can also be found in the form of a series expansion in a. To this end, we 
rewrite (1271) as 



x n - a} j G n , m x m = ay S nfi , 



m=0 



where 



G_ h2n+2m+2 
71.. m. ^ T r)..m.'' / i 



y = a y, 



(36) 

(37) 
(38) 
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or in operator form (see Appendix for notation) 

i" — aG\ x = ay. 
If ||aCr|| < 1, solution of this equation can be represented as a series 

r 1 °° 

'i-a&Y y = ^a n+l G n y. 

71=0 

From here we obtain expansion for Xq 

00 00 
x = ay + a 3 2J G 0tk G kj0 y + a 4 2J G , n G n , k G kfi y 



x — a 1 — 



(39) 



(40) 



fc=i 



00 00 



n,fc=l 



+ « 5 S ^ G o.« G ™. feG ^ G ^^ + °( a6 )- ( 41 ) 

m,n=l k =0 

Evaluating the second term using (1291) and (1371) and comparing this expression with (I34|) we 
derive an alternative expression for Mh 

M h = I + a 2 lJ2( 2k + l)\S k+ i\ 2 h 4k+4 + a 3 J2 G , n G n , k G k , 



k=l 



n,k=l 



00 00 



+ a 4 Go, n G n , k G ktrn G mi o + Of 



a" 



(42) 



m,n=l fc=0 

If all lattice sums (!24j) are real then expression for Af^ is reduced to 



M h = I + a 2 lJ2( 2k + l ) S l+i h 

k=l 

, ^ (2n + 2k + 1) 
+ « 2-) 



4fc+4 



n,k- 



' (2n)!(2fc)! 



5'n+l'S'rt+fc+l'S'/ 



fc+1 



1 
-1 



^4n+4/c+6 



00 00 



1 4r ^ ^ (2rc + 2fc + l)!(2fc + 2m + l). c , 

/ 4 / j m__\\im_\\im_ , iM/n.__M ^n+l >->n+k+l >Jm+k+l <~>m+l '* 



4m+4n+4A;+8 



+ 0(a b ) 



m,n=l k=0 
5> 



(2ra)!(2fc)!(2fc + l)!(2m)! 



(43) 



III. AVERAGE FIELD 



Let us find the average electric field (E) in the parallelogram ABCD in Figured] of area 
S = 4uilmuj2. 

(E) = ^ I EdS = \ I E ^dS + ^ I E ex dS, (44) 

J S J Si n J Sex 



where Si n = ^a 2 is the area of inclusion's cross-section and S ex = S — Si n . From (TIB"]) and 
OlSp we find the average field in the inclusion in Cartesian coordinates 

1 



(E in ) 



S 



E„, ,W - ^ [-A , B V = _^4_ [Co , D )\ 



^ i ri 



S 



\P%n ®ex)S 



(45) 



To calculate the average field outside the inclusion, observe that by Green's theorem 

T 



E ex dS 





du ex 


du ex 




dx 


dy _ 



dS 



<b \-u ex dy, u ex dx] T - <b [-u ex dy, u ex dx] T . 

J ABCD JL 



(46) 



Using the quasiperiodicity property ©-(JH]) of the C-function, one can simplify the integrals 



pB pD pC 

u ex dy = u ex dy+ u ex dy = Re / (-$ ex (z) + $ ex (z + 2wi)) dy 

ABCD J A JC JD 

Re (2u l (-E x + iE y ) + 2Ea 2 {C + iDq)^) Im (2w 2 ). 



p pB pC pD pA 

d> u ex dx = u ex dx + J u ex dx + I u ex dx + I u ex dx 

J ABCD J A JB JC JD 



= Re f ($ ex (z + 2wi) - <$>ex(z)) dx - Re f ($ ex (z + 2uj 2 ) - dx 

JD J D 

= Re ((-E x + iE y )2u l + 2Ea 2 {C + iDq)^) Re (2w 2 ) 
- Re ((-E x + iE y )2u 2 + 2Ea 2 (C + iD )rj 2 ) 2u x . 



(47) 



(48) 



f u ex dy = Re (h § ex (z) dy = Re f $ in (z) dy 

JL JL JL 

°° p2tt 

= Re Ea 2 V {A n + iB n ) I [cos(2n + 1)0 + i sin(2n + 1)0} cos tfdtf = 7ra 2 EA . (49) 

n=0 ^° 

Similarly, 



it PT dx = Re ® $pt(z) dx = Re ® ^^fz) dx 



OO „27T 

-Re £a 2 V (A n + iB n ) / [cos(2n + 1)0 + i sin(2n + 1)0] sin = na 2 EB . (50) 

n=0 ^° 



Thus, the average electric field (E) has the following components: 



Ea? Ea? 
\E) X = E X Re [(Co + iD ) m ] = E x (C Rer/i - D^lm^) , (51) 

„ _ 2 Reu 2 Re [(C + iD Q ) r} X \ - wiRe [(C + i-D ) %] 

\tj)y = Ely + Eld 



Ey + Ea' 



C Im?7iIma;2 — -Do (f — Rer^Imco^) 



(52) 



where we have used the Legendre identity (jUJ). 

Using expression (|34|) for the coefficients Cq and D we rewrite ( }5"T|) -( l5"2"|) in matrix form 



(E) 



[E) x 

[E)y 



, 2a 2 a „ r 
I - — *M / ( 



E., 
E„ 



(53) 



where 



Re r/i Im 2o>2 — Im 771 Im 2co>2 

— Im r/ilm 2cj 2 n — Re r^Im 2a>2 
Calculation of the average electric field in the inclusion and in the medium then gives 



(54) 



\Einj 



(E ex ) 



2ixa 2 <j e 



{(T in + 0~ex)S 



M h 



E x 
E„ 



27ra 2 a ex E 



&ex)S 



2%a 2 a ex _ , 2a 2 a 

M h - — — #M ft 



Cira + &ex)S 



s 



E x 
E„ 



£ - (£ m > - — - *:ro- 



(55) 



(56) 



IV. CALCULATION OF THE EFFECTIVE CONDUCTIVITY 



Effective conductivity of an arbitrary periodic array of inclusions is a tensor cr*, 



® x ® xy 



a xy °y 



It relates the average current density (j) and the average electric field (E) 



(57) 



(j)=<T*(E) or 



j)y = Ky(E)x + <(E) y . 



(58) 
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Observe that 



^ = ~S I •* dS = ~S I j mdS + \ I J ' ex dS 

J S J Sin J Spcr. 



(J; 



s 



cr e . 



/ E in dS + -^ / E ex dS 



S 



'in 



(59) 



To determine the effective conductivity, we apply first a unit electric field in the x-direction: 
E = [E x , E y ] T = [1, 0] T . The corresponding averaged vectors are denoted by the superscript 
(1,0). 



(i>. 



tiTv 



(1,0) 
x 

(1,0) 



a a 

x xy 



(E)x 



(1,0) 

x 

(1,0) 



(60) 



Similarly, we apply then the electric field in the y-direction E = [E x ,E y ] T = [0, 1] T and 
compute the averaged field and current density 



Equation (|60|) - (!6T|) can be written as one matrix equation 



a a 

x xy 



a xy a y 



E)f l) 



(61) 



^ °xy 



a xy a y 



(E)i lfi) (E}x 
(E)y m (E) 



(0,1) 



(0,1) 

y 



(0,1) 



(62) 



Thus, 



< 













(3)y lfi) <i>? W 



(E) x lfi) (E) x 
{ E) y ^ (E) y ^ 



(o,i) 



(63) 



Calculation of the electric field and the current density matrices gives 

do 

IV 



2aa 2 



U)S*> (3)y°> l) 



&exl + 



2aa 2 a P 



S 



(irl - *) M h 



(64) 



(65) 



Substituting these expressions in (163]) we obtain the effective conductivity tensor 

a* = a ex (I + n5M h (I-5*M h y 1 ), (66) 
11 



where 

* = ^ (67) 

is proportional to the fractional part of the inclusions. If ||5^Afh|| < 1 then (J — S^fMhf 1 
can be expanded in a convergent series 

oo 

a* = a ex (l + tt6 (*M,) n ) . (68) 

n=0 

Thus, for a lattice with periods 2u\ and 2u; 2 (see Figure [1]) the expansion of the conductivity 
tensor in terms of volume fraction / of inclusions has the form 



In 2 / 2 ,, , kl , 8a 3 _f 3 

r2 



a* x = a ex l + 2a/ + ^ReC(wi)Im2o; 2 + f- |C(wi) | W2wa + 0(f) , (69) 



< = ImC(a; 1 )Im2a; 2 (l + 2a/ + 0(f)) , (70) 

/ 4a 2 f 2 

a* y = a e J 1 + 2a/ + (ix - Re C(u>i)Im 2u 2 ) 

+ (> 2 - 27rReC(a;i)Im2a; 2 + |CM| 2 Im 2 2^ 2 ) + 0(f)^j , (71) 

where a = — — In what follows we use series (1681) for analytic expression of the 
effective conductivity tensor for specific lattices. 



V. EFFECTIVE CONDUCTIVITIES OF SOME LATTICES 
A. Square lattice 

For the square lattice we put 2u>i = £, 2u 2 = il (see Figure [2]). Then one can find^ that 

m 

and from (I54p and (IBTj) we obtain 

All lattice sums f )24|) are real with the only non-zero being S 2 k, k = 1,2,.... Substituting 
these parameters into (1681) we obtain the effective conductivity tensor of the square lattice 

oo 

<T* = a ex {l + 2j2(ufM h ) n+1 }, (74) 

n=0 



7T 

2t 



7T„ V* = -2i> 



tv 2aa 
2 1 ' 5 = ~P- 



(72) 



(73) 
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} 


r 


» 4 

^ r 


V i 

2co 2 










-0 (- 


3 ^1 





FIG. 2. Square lattice of inclusions of radii a with periods 2cji = ^ and 2oj2 = 



where / = — — is the volume fraction of the inclusions. Calculation of matrix in ( 135]) 
gives 



M fc = (1 + 3a 2 S*h 8 + (9a 4 5 2 4 + 7a 2 S 2 ) h 1& ) I + 0(/i 24 ). 
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(75) 



Effective conductivity tensor of the square lattice is isotropic, <x* = cr*J, and for er* we 
obtain from (1741) 



a* = a ex \ 1 + 2«/ + 2« 2 / 2 + 2a 3 / 3 + 2a 4 / 4 



. „ 9a 5 S 4 + 157r 4 a 7 S 2 + 7a 3 S 2 \ „ 
+ 2 a 9 + - s— 2 1 / 9 



7T 8 



j a io + 27a 6 g 4 + 187r 4 a 8 ff 2 + 14a 4 £f j /10 



, , ( 5WS% + 2l7r 4 a 9 5 2 + 21a 5 Sj \ 

^8 



+ 2 | q11 + i ^ 4 j 7 a 



/ 12 90a 8 g 4 + 247r 4 a 1 % 2 + 28o^ 12 

+2 r + ^ r 

+ 0(/ 13 ). (76) 

E, 1 r-v 1 
- = 3.15121, Si= y -r- = 4.25577 correct to five decimal 

(m + vnr ^— ' (m + zra) 8 

n,m n,m ' 

places. Expression ( J76i) is in agreement with known results 1 - 12 ". 



B. Regular triangular lattice 

The effective conductivity tensor is also isotropic in the case of a regular triangular lattice 

(see Figure [3]). Similar to the previous case we put 2u± = £, 2u2 = £e n ^ 3 . Then we find 9 
that 

" 1 = vr *■ = us' (77) 

and result, 

7r . 2aa 2 . . 

*=^, S= — (78) 

All lattice sums (124"]) are real with the only non-zero being S 3 k, k = 1,2, Substitut- 
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FIG. 3. Regular triangular lattice of inclusions of radii a with periods 2ui\ = £ and 2to>2 = £e ln ^ 3 . 

ing these parameters into ( )68|) we obtain the effective conductivity tensor of the regular 
triangular lattice 

oo 

<x* = a ex [l+ 2J2 (afM h ) n+1 Y (79) 

n=0 

itcl 2 

where / = — — is the fractional part of the inclusions. Matrix Mh found from f )35|) is 

M h (h) = (1 + 5a 2 S 2 h 12 ) I + O (h 24 ) . (80) 

The effective conductivity tensor of the regular triangular lattice is isotropic, er* = a* I, 
and for a* we obtain from f[T9"j) 

a* = a ex I 1 + 2a f + 2a 2 f 2 + 2a 3 f + 2a 4 f 4 + 2a 5 / 5 



7 , 135a 3 S? 
32tt 6 



+ 2a e r+(2« 7 + "7r 3 \f 



. . 135a 4 S 2 \ „ 8 / q 405a 5 S 2 \ , 
„, , 135a 6 Sf\ , 10 , / , 675a 7 Sf\ , n 
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2co 2 




— t 
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1 t 







FIG. 4. Rectangular lattice of inclusions of radii a with periods 2ui = 21 and 2ui2 

The latter formula agrees with calculation in Ref. O. Here S3 = 

correct to five decimal places. 



n,m m + ne 3 



5.86303 



C. Rectangular lattice 

Consider a rectangular lattice generated by the vectors 2u>\ = 2£, 2u 2 = it. We compute 
the lattice sums 



(2u)\m + 2u 2 n) A ~ (2m + in) 4 



2.16646, 



S 3 = V' jr = -2.03111. 

^ (2m + zn) 6 



Next, we compute the invariants g 2 and g 3 : 

g 2 = 60 S 2 t~ 4 = 129.988T 4 , 
# 3 = 140 S 3 t~ 6 = -284.355r 6 . 
Then we compute r]i from (JTj) using ffTB"]) 

rya = 92, £3) = r J C (1; 129.988, -284.355) 
= -0.1480OT 1 



(82) 
(83) 



(84) 
(85) 



(86) 
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and matrices ^ and Mh 



-0.14800 



3.28959 



(87) 



M h = (1 + 3a 2 Sih 8 + 5a 2 S 2 h 12 ) I + 30a 3 S 2 S 3 /i 14 



1 



-1 



+ 0{h 
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From (1681) we obtain that the effective conductivity tensor cr* is diagonal 



°* y 



with components 



where / is the volume fraction of the cylinders / = \-kY? . 



(89) 



a * x = a ex (1 + 2af - 0.188439a 2 / 2 + 0.0177547a 3 / 3 - 0.00167284a 4 / 4 + O (/ 5 )) , (90) 
a* = a ex (1 + 2a/ + 4.18844a 2 / 2 + 8.77150a 3 / 3 + 18.3694a 4 / 4 + O (f)) , (91) 



D. Anisotropic lattice 

Here we show how to find the effective conductivity tensor for an arbitrary lattice. Con- 
sider the case when the lattice is created by the vectors 2ui = 2£, 2uj 2 = £e m ^ 3 (see Figure 
[5]). Then we calculate the lattice sums 

d>4 , J 



S2 = £' 



n,m 



S3 = E' 



(2ouim + 2oj 2 n 



I 6 

2m + e^n 



/'./" i 2m + e^n 



-1.08720+ 1.88309i, 



2.01542, 



(92) 



(93) 



and the invariants 



g 2 = 60S 2 r 4 = (-65.2321 + 112.985*) £~\ 
g 3 = U0S 3 r 6 = 282.158 r 6 . 



(94) 
(95) 
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FIG. 5. Generic lattice of inclusions of radii a with periods 2lj% = 2£ and 2co>2 = £e m ^ 3 . 



In order to find constant 7/1 we use the homogeneity property of the (-function ([13 
Vl = C(ux,g 2 , g 3 ) = (1, -65.2321 + 112.985z, 282.158) 
= (1.07651 - 1.277032)r 4 . 
Now we evaluate matrix \I/ in ( 154)) 



(96) 



0.93228 1.10594 
1.10594 2.20931 



(97) 



Re r\\ Im 2o;2 — Im r)\ Im 2o>2 
-Im T^Im 2a>2 n — Re r^Im 2u 2 

and using the expansion of Mh{K) in ( 1351) we compute from ( 1681) the effective conductivity 
tensor 



a* a* 

" x xy 



a xy °y 



(98) 



where 



a* = a ex (1 + 2af + 1.18702a 2 / 2 + 1.69592a 3 / 3 + 2.98936a 4 / 4 + O (f 5 )) , (99) 
a* xy = a ex (1.40813a 2 / 2 + 2.81625a 3 / 3 + 5.15506a 4 / 4 + O (f 5 )) , (100) 
a* = a ex (1 + 2a/ + 2.81298a 2 / 2 + 4.94784a 3 / 3 + 8.94191a 4 / 4 + O (f 5 )) . (101) 



APPENDIX 
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Here we study properties of system (|27|) using the approach similar to that in Ref. 
We seek a solution of ( 127]) in the space of bounded sequences /^(M 2 ) whose elements are 
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two-dimensional vectors 

where x n = [xn , xi 2 '] T . The norm of an element of this space is given by 

||cc|| = sup max {| 2^|, |} ■ 

n=0,l,2,... 

We introduce a linear operator @(h) by 

oo 

^(h)x) n = Y,G n , m x m h 2n+2m+2 , n = 0,1,2,.... (Al) 

m=0 

Then (j2"Tj) can be written in operator form 

x-&(h)x = y. (A2) 
Properties of operator &(h) and equation (1A2I) are summarized in the following 

Theorem 1. For each ^ h ^ ~ vis a bounded operator in /^(R 2 ). If0^.h<h then 

@(h) is compact and can be represented by a convergent series 

oo 

^(/i) = ^/i 2m G (m) , (A3) 

771=1 

where G^ are finite- dimensional operators of order 2m. 

Proof. Let us estimate the norm of @(h): 

\\&(h)x\\ II V°° G x h 2n + 2m + 2 \\ 

\\&(h)\\ = SUP " " = SUp SUp 11 ^m=0^n,mX m tl || 



||as||<l \\ x \\ IklKl n H^ll 

oo 

< sup ]T ||G n , m || h 2n+2m+2 <: \a\ sup (\S*\ + |5j|) 

n n n 

m=U 

X Sim V ( 2n + 2m + X ) ! fc2n+2m+2 _ I I rf /^2«+2 ^n+l °^ +2n+1 

X S ? 2. (2 TO )! (2n + 1)! " H 6 T (2n + 1)! rf^ 1 2- 

m=0 v ' v 7 m=0 

^n+2 rf 2n+l 

= |a| b sup 



n (2n + l)!rf/i 2n+1 V 1 + ^ 



, \ 2n+2 / , \ 2n+2 N 

h \ in 



a\ S sup + 



l-h \l + h 



h V ( 



\«\s[(— h + l <0 °' (A4) 
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where S = sup n (\S^\ + |<S^|). Therefore @(h) is a bounded operator for ^ h ^ |. From 
flAil) it also follows that if ^ h < ~ than ^(/i) maps a bounded sequence into the space 
c (M 2 ) of sequences converging to zero and hence it is compact. 

Expansion ( 1A3I) follows formally from the definition (1A1I) of operator &(h), where 2m- 
dimensional operators are defined by 

{G m - k -i ik x k , < A; < m - 1, 
(A5) 
0, k ^ m. 

To show convergence of the series (1A3j) we observe that 

< M § P"'- 1 )' (A6) 
1 1 ml (m- 1)! v ; 

Therefore series ( 1A3I) is dominated by a convergent for ^ h < | series 

^i" ^m!(m-l)! (l + vT=4/?) 1 ; 

For values of a, S, and h such that < 1 in (1A4j) the fixed point theorem for 

contraction operators on Banach spaces ensures the following properties of the solution of 

m 



(a) Equation ( 1Z7)) has a unique solution aj G c (K 2 ). 

(b) The truncated solution of (|27|) converges exponentially to a? . 

(c) The solution of ( 1271) can be represented as a convergent power series in h. 

□ 
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